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Abstract. Numerical relativity has made big strides over the 
last decade. A number of problems that have plagued the 
field for years have now been mostly solved. This progress 
has transformed numerical relativity into a powerful tool to 
explore fundamental problems in physics and astrophysics, and 
I present here three representative examples. These "three 
little pieces" reflect a personal choice and describe work that I 
am particularly familiar with. However, many more examples 
could be made. 



1. Introduction 

Numerical relativity has hardly seen better times before. Over the last 
few years, in fact, a truly remarkable development has shaken the field. 
Starting from the first simulations showing that black-hole binaries could 
be evolved for a few orbits [111) \35\ [23] , or that black- hole formation could 
be followed stably using simple gauges and without excision [22 , new 
results, some awaited for decades, have been obtained steadily. As a direct 
consequence of this "Renaissance", it is now possible to simulate binary 
black holes [42 and binary neutron stars accurately for dozens of 
orbits, from the weak-field inspiral, down to the final black-hole ringdown 
(see also [38 l 1107] for recent reviews). 

There are several reasons behind this rapid progress. These include 
the use of more advanced and accurate numerical techniques |116[ 1117] , 
the availability of larger computational facilities, but also the development 
of formulations of the Einstein equations and gauges that are particularly 
well-suited for numerical evolutions [95l l3l [T35l 12^ 11 [TT2l l5l [27l IM^^ 
[6]. The phase transition that has taken during this year has radically 



changed numerical relativity, freeing it from the corner of idealised 
investigations. Most importantly, it has transformed numerical relativity 
into a research area where long-standing problems can found a quantitative 
and accurate solution, and into a tool by means of which it is possible to 
explore fundamental aspects of physics and astrophysics. 

Numerous examples could be given to testify this transformation, 
although I will report here only those that I am particular familiar 
with. More specifically, in what follows I will discuss: (i) how numerical 
simulations of magnetised neutron stars provide convincing evidence 
that this process leads to the conditions that are expected behind the 
phenomenology of short gamma ray burst; (ii) how numerical simulations 
of the head-on collision of selfgravitating fiuids boosted at relativistic 
speeds can be used to understand the conditions leading to the formation 
of a black hole and provide a dynamical version of the hoop conjecture; 
(in) how the study of the local properties of apparent horizons can be 
used to explain bizarre behaviours in binary black-hole simulations and 
can be effectively correlated with a portion of the spacetime infinitely far 
away: This selection is by no means comprehensive, but rather a 

very personal one, and I apologise in advance for not discussing all the 
excellent work that cannot find space in this contribution. 

2. First piece: Prom neutrons star to gamma-ray bursts 

The numerical investigation of the inspiral and merger of binary neutron 
stars in full general relativity has seen enormous progress made in recent 
years. Crucial improvements in the formulation of the equations and 
numerical methods, along with increased computational resources, have 
extended the scope of early simulations. These developments have made 
it possible to compute the full evolution, from large binary-separations up 
to black-hole formation, without and with magnetic fields |137[ [T9| [H [86l 
[SHI |60] , and with idealised or realistic equations-of-state [121" TT . This 
tremendous advancement is also providing information about the entire 
gravitational waveform, from the early inspiral up to the ringing of the 
black hole (see, e.g., [181 [Ml l2Q])- Advanced interferometric detectors 
starting from 2014 are expected to observe these sources at a rate of 
~ 40 — 400 events per year [1 . 

These simulations also probe whether the end-product of mergers can 
serve as the "central engine" of short gamma-ray bursts (SGRBs) |103[ 
[53l [97] . The prevalent scenario invoked to explain SGRBs involves the 
coalescence of a binary system of compact objects, e.g., a black hole and a 
neutron star or two neutron stars [129, 128, 96, 84 . After the coalescence, 
the merged object is expected to collapse to a black hole surrounded by an 
accretion torus. An essential ingredient in this scenario is the formation 



of a central engine, which is required to launch a relativist ic outflow with 
an energy of ~ 10^^ — 10^^ erg on a timescale of ~ 0.1 — 1 s [96^ 84 . 

The qualitative scenario described above is generally supported by the 
association of SGRBs with old stellar populations, distinct from the young 
massive star associations for long GRBs [56^, 115" . It is also supported to a 
good extent by fully general-relativistic simulations, which show that the 
formation of a torus of mass Mtor ^ 0.4 Mq around a black hole with spin 
J/M^ 0.7 — 0.8, is inevitable |121j . In addition, recent simulations have 
also provided the first evidence that the merger of a binary of modestly 
magnetised neutron stars naturally forms many of the conditions needed 
to produce a jet of ultrastrong magnetic field, with properties that are 
broadly consistent with SORB observations. This missing link between the 
astrophysical phenomenology of GRBs and the theoretical expectations is 
a genuine example of the new potential of numerical relativity and I will 
discuss it in detail belo^vfl 

2.1. The numerical setup 

It is not useful to discuss here in detail the numerical setup and the 
technical details of the numerical codes used in these calculations. These 
details can be found in Refs. [57l |60], while a description of the physical 
initial data was presented in Ref. [122^. It is sufficient to recall here 
that the evolution of the spacetime is obtained using a three-dimensional 
finite-differencing code providing the solution of a conformal traceless 
formulation of the Einstein equations [109^ {i.e., the CCATIE code). The 
equations of general-relativistic magnetohydrodynamics (GRMHD) in the 
ideal-MHD limit are instead solved using a code code [20l [211 IHZ] which 
adopts a flux-conservative formulation of the equations as presented in [11] 
and high-resolution shock-capturing schemes {i.e., the Whisky code). In 
order to guarantee the divergence-free character of the MHD equations 
the flux- CD approach described in jl44) was employed, although with the 
difference that the vector potential is used as evolution variable rather 
than the magnetic field. Both the Einstein and the GRMHD equations 
are solved using the vertex-centred adaptive mesh-refinement (AMR) 
approach provided by the Carpet driver [131 . In essence, the highest- 
resolution refinement level is centred around the peak in the rest-mass 
density of each star and in moving the "boxes" following the position 
of this maximum as the stars orbit. The boxes are evolved as a single 
refinement level when they overlap. The calculations were carried out 
using six levels of mesh refinement with the finest level having a resolution 
of A = 0.1500 M© 221m. 

^ Much of what follows is taken from the discussion presented in Ref. [122| . 



From a physical point of view, the initial data is given by a configuration 
that could represent the properties of a neutron star-binary a few orbits 
before their coalescence. More specifically, we simulate two equal-mass 
neutron stars, each with a gravitational mass of 1.5 Mq {i.e., sufficiently 
large to produce a black hole soon after the merger) , an equatorial radius of 
13.6 km, and on a circular orbit with initial separation of o:^ 45 km between 
the centres (all lengthscales are coordinate scales). Confined in each star 
is a poloidal magnetic field with a maximum strength of 10^^ G. At this 
separation, the binary loses energy and angular momentum via emission of 
gravitational waves, thus rapidly proceeding on tighter orbits as it evolves. 

2.2. The basic dynamics 

After about 8 ms (~ 3 orbits) the two neutron stars merge forming a 
hypermassive neutron star (HMNS), namely, a rapidly and differentially- 
rotating neutron star, whose mass, 3.0 M©, is above the maximum mass, 
2.1 M0, allowed with uniform rotation by an ideal-ffuid equation of state 
(EOS)r] p = {T — l)pe, where p is the baryonic density, e the specific 
internal energy, and F = 2 with an adiabatic index of 2. Being metastable, 
a HMNS can exist as long as it is able to resist against collapse via a 
suitable redistribution of angular momentum [e.g., deforming into a "bar" 
shape |137[ [19]]. or through the increased pressure-support coming from 
the large temperature-increase produced by the merger. However, because 
the HMNS is also losing angular momentum through gravitational waves, 
its lifetime is limited to a few ms, after which it collapses to a black hole 
with mass M = 2.91 Mq and spin J/M^ = 0.81, surrounded by a hot and 
dense torus with mass Mtor = 0.063 Mq [60 . 

These stages of the evolution can be seen in Fig. [l] which shows 
snapshots of the density colour-coded between 10^ and lO^^gr/cm^, and 
of the magnetic field lines (green on the equatorial plane and white outside 
the torus). Soon after the black hole formation the torus reaches a 
quasi-stationary regime, during which the density has maximum values of 
~ 10^^ g/cm^, while the accretion rate settles to M ~ 0.2 Mq/s. Using the 
measured values of the torus mass and of the accretion rate, and assuming 
the latter will not change significantly, such a regime could last for taccr — 
Mtor/M ~ 0.3 s, after which the torus is fully accreted; furthermore, if 
the two neutron stars have unequal masses, tidal tails are produced which 
provide additional late-time accretion [121 . This accretion timescale is 
close to the typical observed SGRB durations (HTJ [96]. It is also long 
enough for the neutrinos produced in the torus to escape and annihilate in 
its neighbourhood; estimates of the associated energy deposition rate range 

^ The use of a simplified EOS does not infiuence particularly the results besides 
determining the precise time when the HMNS collapses to a black hole. 



Figure 1. Snapshots at representative times of the evolution of the binary and of the 
formation of a large-scale ordered magnetic field. Shown with a colour-code map is the 
density, over which the magnetic-field lines are superposed. The panels in the upper 
row refer to the binary during the merger (t = 7.4 ms) and before the collapse to black 
hole {t = 13.8ms), while those in the lower row to the evolution after the formation of 
the black hole {t = 15.26ms, t = 26.5 ms). Green lines sample the magnetic field in the 
torus and on the equatorial plane, while white lines show the magnetic field outside the 
torus and near the black hole spin axis. The inner/outer part of the torus has a size of 
^ 90/170 km, while the horizon has a diameter of '2^ 9 km. 



from ~ 10^^ erg/s |48] to ~ 10^^ erg/s [134j . thus leading to a total energy 
deposition between a few 10^^ erg and a few 10^^ erg over a fraction of a 
second. This energy would be sufficient to launch a relativistic fireball, but 
because radiative losses are accounted yet, the large reservoir of thermal 
energy in the torus cannot be extracted in these simulations. 

The gravitational wave signal of the whole process is shown in the 
top part of the left panel in Fig. [2j while the bottom part exhibits the 
evolution of the MHD luminosity, Lmhd , as computed from the integrated 
Poynting flux (solid line) and of the corresponding energy, £^mhd, (dashed 
line). Clearly, the MHD emission starts only at the time of merger and 
increases exponentially after black- hole formation, when the gravitational 




Figure 2. Left panel: gravitational wave signal shown through the i = 2,m = 2 mode 
of the + polarization, (/i+)22, (top part) and of the MHD luminosity, Lmhd? (bottom 
part) as computed from the integrated Poynting flux and shown with a solid line. The 
corresponding energy, ^^mhd, is shown with a dashed line. The dotted and dashed 
vertical lines show the times of merger (as deduced from the flrst peak in the evolution 
of the gravitational wave amplitude) and black-hole formation, respectively. Right 
panel: Evolution of the maximum of the magnetic fleld in its poloidal (red solid line) 
and toroidal (blue dashed line) components. The bottom panel shows the maximum 
local fluid energy indicating that an unbound outflow {i.e.^ -E'loc > 1) develops and is 
sustained after black-hole formation. 



wave signal essentially shuts off. Assuming that the quasi-stationary MHD 
luminosity is ~ 4 x lO^^erg/s, the total MHD energy released during the 
lifetime of the torus is ~ 1.2 x 10^^ erg, which, if spread over an opening 
half-angle of ~ 30° (see discussion below), suggests a lower limit to the 
isotropic equivalent energy in the outflow of ~ 9 x 10^^ erg. While this is at 
the low end of the observed distribution of gamma-ray energies for SGRBs, 
larger MHD luminosities are expected either through the additional growth 
of the magnetic field via the winding of the field lines in the differentially- 
rotating disk (the simulation covers only one tenth of taccr), or when 
magnetic reconnection (which cannot take place within an ideal-MHD 
approach), is also accounted for [which may also increase the gamma-ray 
efficiency, e.g., [92 J. 

The last two panels of Fig. [l] offer views of the accreting torus after the 
black-hole formation. Although the matter dynamics is quasi-stationary, 
the last two panels clearly show that the magnetic-field is not and instead 
evolves significantly. It is only when the system is followed well after the 
formation of a black hole, that MHD instabilities develop and generate 
the central, low-density, poloidal-field funnel. This regime, which was 



not accessible to previous simulations [ 114[ |9l [86], is essential for the 
jet formation [79j. Because the strongly magnetised matter in the 
torus is highly conductive, it shears the magnetic-field lines via differential 
rotation. A measurement of the angular- velocity in the torus indicates that 
it is essentially Keplerian and thus unstable to the magneto-rotational 
instability [25 , which develops :^ 5 ms after black- hole formation and 
amplifies exponentially both the poloidal and the toroidal magnetic fields; 
the e- folding time of the instability is 2.5 ms and in good agreement 
with the one expected in the outer parts of the torus Because of 

this exponential growth, the final value of the magnetic field is largely 
insensitive to the initial strength and thus a robust feature of the dynamics 
(see also [138 for a similar behaviour recently computed in a HMNS) 

A quantitative view of the magnetic-field growth is shown in the right 
panel of Fig. [2j which shows the evolution of the maximum values in the 
poloidal and toroidal components. Note that the latter is negligibly small 
before the merger, reaches equipartition with the poloidal field as a result 
of a Kelvin-Helmholtz instability triggered by the shearing of the stellar 
surfaces at merger [114, 59 , and finally grows to 10^^ G by the end 
of the simulation. At later times {t > 22 ms) , when the instability is 
suppressed, the further growth of the field is due to the shearing of the 
field lines and it increases only as a power-law with exponent 3.5 (4.5) for 
the poloidal (toroidal) component. Although the magnetic-field growth 
essentially stalls after t 35 ms, further slower growths are possible |100] , 
yielding correspondingly larger Poynting fiuxes. Indeed, when the ratio 
between the magnetic fiux across the horizon and the mass accretion 
rate becomes sufficiently large, a Blandford-Znajek mechanism [31 may 
be ignited [78]; such conditions are not met over the timescale of the 
simulations, but could develop over longer timescales. Also shown in the 
right panel of Fig. |2] is the maximum local ffuid energy, highlighting that 
an unbound outflow {i.e., E\oc > 1) develops after black-hole formation 
along the outer walls of the torus and persists for the whole duration of 
the simulation. 

Finally, Fig. [3| provides a summary of the magnetic-field dynamics. It 
shows the magnetic field in the HMNS formed after the merger and its 
structure and dynamics after the collapse to black hole. In particular, in 
the last three panels it shows the magnetic-field structure inside the torus 
and on the equatorial plane (green), and outside the torus and near the 
axis (white). It is apparent that the highly turbulent magnetic field in the 
HMNS {t = 13.8 ms) changes systematically as the black hole is produced 
{t = 15.26 ms), leading to the formation of a toroidal magnetic field in the 
torus {t = 21.2 msj^ As the MRI sets in, the magnetic field is not only 

^ Turbulence in relativistic flows is an extremely challenging problem that is also 



Figure 3. Magnetic-field structure in the HMNS (first panel) and after the collapse 
to black hole (last three panels). Green refers to magnetic-field lines inside the torus and 
on the equatorial plane, while white refers to magnetic-field lines outside the torus and 
near the axis. The highly turbulent, predominantly poloidal magnetic-field structure 
in the HMNS (t = 13.8 ms) changes systematically as the black hole is produced 
(t = 15.26 ms), leading to the formation of a predominantly toroidal magnetic field 
in the torus (t = 21.2ms). All panels have the same linear scale, with the horizon's 
diameter being of ^ 9 km. 



amplified, but also organises itself into a dual structure, which is mostly 
toroidal in the accretion torus with 5tor ~ 2 x 10^^ G, but predominantly 
poloidal and jet-like along the black-hole spin axis, with B^oi ~ 8 x 10^^ G 
{t = 26.5 ms|^ Note that the generation of an ordered large-scale field is 
far from trivial and a nonlinear dynamo may explain why the MRI brings 
a magnetic field self-organization, as it has been also suggested in case of 
MRI-mediated growth of the magnetic field in the conditions met in the 

essentially unexplored. Also in this case, the first relativistic simulations have been 
performed only recently [1181 1150) . 

^ A similar magnetic-field configuration has been recently reproduced also when 
simulating the merger of a magnetised neutron star onto a black hole ^54j. 



collapse of massive stellar cores |85 l llQQj . However, the jet-like structure 
produced in the simulation is not yet the highly collimated ultrarelativistic 
outflow expected in SGRBs (see also below). 

The magnetic hollow jet-like structure has an opening half-angle of 
~ 30°, which sets an upper limit for the opening half-angle of any 
potential outflow, either produced by neutrino energy deposition [7 or by 
electromagnetic processes [79]. In these simulations most of the outflow 
develops along the edges of the jet-like structure, via a turbulent layer 
of electromagnetic driven matter, which shields the central funnel from 
excessive baryonic pollution. It is reasonable to expect that such a layer is 
crucial to set the opening angle of any ultrarelativistic jet, to shape both 
the radial and transverse structure of the jet, as well as to determine its 
stability properties. The Lorentz factors of the outflow measured in these 
simulations are not very high (7 < 4), but can potentially be amplified 
by several orders of magnitude in the inner baryon-poor regions through 
special-relativistic effects [8 , the variability of the flow [64], or when 
resistive-MHD effects are taken into account [50 . Such accelerations will 
be produced as a more realistic and general-relativistic treatment of the 
radiative losses will become computationally affordable. 

2.3. Comparison with observations 

Below I briefly discuss how the results presented above broadly match the 
properties of the central engine as deduced from the observations. 

Duration: The observed duration of the prompt gamma-ray emission 
GRBs is energy dependent and is usually determined through T^^, the time 
over which x% of the total counts are observed, between the (100 — x)/2 
and (100 + x)/2 percentiles. The most common intervals used are T90 
(or T50), initially defined ^ between 20 keV and 2 MeV. The ORB 
duration distribution is bimodal [81], where the durations of SGRBs 
(approximately 25% of GRBs) are well-fit by a fairly wide log-normal 
distribution centred around T90 ^ 0.8 s with a FWHM of 1.4 dex [96]. 
The typical redshifts of the SGRBs observed with Swift are in the range 
z ~ 0.3 — 1, suggesting a central value of the intrinsic duration distribution 
of ~ (1 + 2:)~^0.8 s ~ 0.5 s, and a comparably wide distribution around 
this value. This is in close agreement with our accretion time of ~ 0.3 s. 

Energy: The isotropic equivalent energy output in the prompt gamma-ray 
emission of SGRBs, E^^i^o^ spans a wide range, from (2.7±1) x 10^^ erg (in 
the observed energy range 15 — 350 keV) for GRB 050509B at a redshift 
of z = 0.225 [33], up to (1.08 ± 0.06) x 10^^ erg (in the observed energy 
range 10 keV - 30 GeV) for GRB 090510 at z = 0.903 [3 . However, the 
most typical values are in the range Ej^i^o ^ 10^^ — 10^^ erg In this 
model, the highly relativistic outflow may be powered either by neutrino- 



anti neutrino annihilation, or by the Blandford-Znajek mechanism. For 
the former one might expect a total energy release between a few 10^^ erg 
and ~ 10^^ erg [101, 29 , into a bipolar relativistic jet of opening half- 
angle 6>jet 8 — 30°, corresponding to a fraction ^ 0.01 — 0.13 of the 
total solid angle, and isotropic equivalent energies, £^j,^^iso, between a few 
10^^ erg and ^ 10^^ erg. For the latter mechanism, instead, and if the 
magnetisation near the event horizon becomes sufficiently high, the jet 
power for these values for the black-hole mass and spin is |83j 

i„,.3.0xl0"(&!)(^^)%.g/s. (1) 

where /rei is the fraction of the total Blandford-Znajek power that is 
channelled into the resulting relativistic jet (and /rei ~ 0.1 might be 
expected for ejecta with asymptotic Lorentz factors above 100). This 
relativistic outflow is launched over a timescale of ~ 0.2 s and corresponds 
to 

£BZ,,..~I.2xI0"(if)(4)"(^^)%,g, (2, 

Comparing the X-ray afterglow luminosity (after 10 or 11 hours) and 
£^7,iso suggests that the efficiency of the prompt gamma-ray emission in 
SGRBs is typically high and similar to that of long GRBs with 
^7,iso (0-1 — 0.9)£^iso, radiating between ~ 10% and ~ 90% of the initial 
energy of the ultrarelativistic outflow. Therefore, this model is able to 
accommodate the observed £^^,iso values. 

Lorentz factor: The Fermi Gamma- Ray Space Telescope has detected 
GeV emission from SGRBs [2 , suggesting typical lower limits of 7min ^ 
10^ - 10^. In particular, 7min^ 1200 was obtained for GRB 090510 [3]. 
However, a more realistic model [63 results in 7min values lower by a 
factor of ~ 3. Therefore, the central engine should be capable of producing 
outflow Lorentz factors of at least a few hundred. The fact that our 
simulation produces a strongly magnetised mildly relativistic outflow at 
angles near ~ 30° from the black-hole spin axis would help shield the 
inner region near the spin axis from excessive baryon loading, and thus 
assist in achieving high asymptotic Lorentz factors at large distance from 
the source, after the outflow in this region is triggered by neutrinos and/or 
the Blandford-Znajek mechanism. 

Jet angular structure: This is poorly constrained by observations (even 
more so than for long GRBs). The only compelling case for a jet break 
in the afterglow light-curve is for GRB 090510 [44], which occurred very 
early on (after ~ 1400 s), and would thus imply an extremely narrow 
jet (^jet ~ 0.2 — 0.4°) and modest true energy output in gamma-rays 
(~ 10^^ erg). If this is indeed a jet break, it might correspond to a 



line of sight near a very narrow and bright core of a jet, which also has 
significantly wider wings. Observers with lines of sight along these wings 
would then see a much dimmer and more typical SGRB [127, 104 ; without 
such wings, however, the observations would suggest a very large intrinsic 
and beaming-corrected event rate per unit volume. In most cases there 
are only lower limits on a possible jet break time [96 , resulting in typical 
limits of /5 > 10~^ or ^jet ^ 8°. This is consistent with the expectation of 
^jet ^ 8 — 30° for the ultrarelativistic ejecta capable of producing a SGRB 
(which would also imply a reasonable SGRBs intrinsic event rate per unit 
volume) . 

2.4- Summary 

The calculations reported above demonstrate that a binary merger of two 
neutron stars inevitably leads to the formation of a relativist ic jet-like 
and ultrastrong magnetic field, which could serve as a central engine for 
SGRBs. Because the magnetic-field growth is exponential, the picture 
emerging from these simulations is rather general and applies equally even 
to mildly magnetised neutron stars. Overall, this first "little piece" of 
numerical relativity removes a significant uncertainty as to whether such 
binary mergers can indeed produce the central engines of SGRBs. While 
the electromagnetic energy release is already broadly compatible with the 
observations, the simulations discussed above lack a proper treatment of 
the energy losses via photons and neutrinos or resistive dissipation, which 
can provide a fundamental contribution to the energy-input necessary 
to launch the fireball and cool the torus jl34[ [48]. This additional 
energy input, whose self-consistent inclusion in general relativity remains 
extremely challenging, may help to launch an ultrarelativistic outfiow very 
early after the black hole forms and complete the picture of the central 
engine of a SGRB. 

3. Second piece: A dynamical hoop conjecture 

The second "little piece" of numerical relativity that I will discuss 
aims at addressing the issue of necessary conditions for the formation 
of a black hole, which still represents one of the most intriguing 
and fascinating predictions of classical general relativity. There is 
abundant astronomical evidence that black holes exist, and a number 
of considerations supporting the idea that under suitable conditions 
gravitational collapse is inevitable [145 . In addition, there is overwhelming 
numerical evidence that black-hole formation does take place in a variety 
of environments [122 . Yet, a rigorous definition of the sufficient conditions 
for black- hole formation is still lacking. Hence, it is not possible to predict 



whether the cohision of two compact objects, either stars or elementary 
particles, will lead to the formation of a black hole. 

The hoop conjecture proposed by Thorne in the '70s, provides some 
reasonable and intuitive guidelines |142j . I recall that the conjecture states 
that a black hole is formed if an amount of "mass-energy" E can be 
compressed to fit within a hoop with radius equal or smaller than the 
corresponding Schwarzschild radius, i.e., if i^hoop ^ Rs = 2GE/c^^ where 
G is gravitational constant and c the speed of light. Even though it can 
be made precise under particular circumstances [133 , the hoop conjecture 
is not meant to be a precise mathematical statement and, in fact, it is 
difficult to predict if the above-mentioned collision will compress matter 
sufficiently to fit within the limiting hoop. Loosely speaking, what is 
difficult is to determine which part of the "kinetic energy" of the system can 
be accounted to fit within the hoop. Since at the collision the conversion 
of kinetic energy into internal energy is a highly nonlinear process, any 
quantitative prediction becomes rapidly inaccurate as the speeds involved 
approach that of light. 

As stated above, the hoop conjecture is purely classical. A quantum- 
mechanical equivalent is not difficult to formulate, although not very 
stringent, as it simply implies that a black hole will be formed at Planck- 
energy scales. The predicting power does not improve significantly when 
considering the conditions of black-hole formation in higher-dimensional 
theories of gravity (see, e.g., [HI I149[ |148^). In these frameworks, the 
energy required for black- hole formation might be significantly smaller [T^ , 
thus providing the possibility of producing them in the Large Hadron 
Collider (LHC) [49 , but no firm conclusion has been reached yet. 

Clearly, although numerical simulations represent a realistic route to 
shed some light on this issue (see, e.g., , |14Q(ll36(ll41j ). even the simplest 
scenario of the collision of two compact objects at ultrarelativistic speeds 
is far from being simple and it is actually very challenging. A ffist step was 
taken by Eardley and Giddings [51 , who have studied the formation of a 
black hole from the head-on collision of two plane-fronted gravitational 
waves with nonzero impact parameter (previous work of D'Eath and 
Payne j45] |46l |47] using different methods had considered a zero impact 
parameter). In all of these analyses each incoming particle is modelled as 
a point particle accompanied by a plane-fronted gravitational shock wave 
corresponding to the Lorentz-contracted longitudinal gravitational field of 
the particle. At the instant of collision the two shock waves pass through 
one another and interact through a nonlinear focusing and shearing. As a 
result of their investigation, a lower bound was set on the cross- section for 
black- hole production, i.e., a > 32.5(G£^/2c'^)^, where E is the centre-of- 
mass (lab) energy. More recently, and in a framework which is closer to the 
one considered here, this problem has been investigated by Choptuik and 



Pretorius [39], who studied the cohision of two classical spherical solitons, 
with a total energy of the system in the lab frame E = 27bmoc^, where mo 
is the "rest-mass" , 7b = 1/ y^l — v^/c^ and vi^ the boost velocity. They 
were then able to show that for collisions with sufficiently high boost, 
^•^v 7b ^ 2.9, a black hole can be formed. 

In what follows I discuss what has been recently reported on the first 
calculations of black- hole production from the collision of two compact, 
selfgravitating, fluid objects boosted at ultrarelativistic speed^ (A similar 
investigation by East and Pretorius [52^ has also appeared at about the 
same time). 

I start by pointing out that there are several important differences 
with the previous investigations in ^| ^0 \^ [5ll 09]. Differently 
from [45l [46l [47l |5T] , in fact, I will consider colliding objects that are not in 
vacuum and are not treated as point particles. Rather, they are relativistic 
stars, which obviously extended and selfgravitating objects, thus with a 
behaviour that is intrinsically different. Also, differently from [39|, these 
objects are not described as scalar fields, but as fluids and thus represent 
a more realistic description of baryonic matter, such as the one employed 
when simulating relativistic heavy-ion collisions |125j . These intrinsic 
differences also make the comparison with the works of [45^, [46l [47l |5l] 
very hard if possible at all. On the other hand, many analogies exist with 
the collision of bosons stars considered in [39 , and that, as I will discuss 
below, can be interpreted within the more general description of black- hole 
production from ultrarelativistic collisions. 

Overall, the most important and distinguishing feature in the collision 
of two selfgravitating stars is that a black hole can be produced even from 
zero initial velocities if the initial masses are large enough; this behaviour 
is clearly absent in all previous results, where instead a critical initial boost 
is necessary [45j [46l [47l EB [39] . In addition, for each value of the effective 
Lorentz factor, (7), a critical initial mass exists, Mc, above which a black 
hole is formed and below which matter, at least in part, selfgravitates. 
More importantly, both Mc follows a simple scaling with (7) , thus allowing 
to extrapolate the results to the masses and energies of modern particle 
accelerators and to deduce that black-hole production is unlikely at LHC 
scales. 

3.1. The numerical setup 

The numerical setup employed in the simulations is the same presented 
in [76], and it uses an axisymmetric code to solve in two spatial dimensions, 
the set of the Einstein and of the relativist ic-hydrodynamic 
equations [75]. The axisymmetry of the spacetime is imposed exploiting 

^ Much of what follows is taken from the discussion presented in Ref. [124| . 



the "cartoon" technique, while the hydrodynamics equations are written 
exphcitly in cyhndrical coordinates. Ah the simulations use an ideal- fluid 
EOS with r = 2. The initial configurations consist of spherical stars, 
constructed as in [76l 1119] after specifying the central density, pc^ where 
the latter also serves as parameter to determine the critical model. The 
stars have an initial separation D and are boosted along the z-direction 
via a Lorentz transformation with boost vi, / c. To limit the initial violation 
in the constraints, D is chosen to be sufficiently large, i.e., D = 240 M©, 
and an optimal superposition of the two isolated-star solutions that will 
be presented in a longer paper. The grid has uniform spacing A = 
0.08(0.06) Mo with extents x/Mq e [0,80] and z/Mq G [0,150(200)], 
where the round brackets refer to the more demanding high-boost cases. 
Refiection boundary conditions are applied on the z = plane, while 
outgoing conditions are used elsewhere. 

3.2. The basic dynamics 

The dynamics of the process is rather simple. As the two stars 
approach each other, the initial boost velocity increases as a result of 
the gravitational attraction, leading to a strong shock as the surfaces of 
the stars collide. In the case of supercritical initial data, i.e., of stars with 
masses above a critical value, Mc, a black hole is promptly produced and 
most of the matter is accreted. Conversely, in the case of suhcritical initial 
data, i.e., of stars with masses below Mc, the product of the collision is 
a hot and extended object with large- amplitude oscillations. Part of the 
stellar matter is unbound and leaves the numerical grid as the product of 
the collision reaches an equilibrium. 

Figure [4] shows snapshots at representative times of the rest-mass 
density, p (top row), of the Lorentz factor, 7 = (1 — v^vij (middle 
row), and of the local fiuid energy, —u^ (bottom row), for two subcritical 
models. The left panel, in particular, refers to a binary boosted at 
v\^/c = 0.3. Note that the stars are strongly compressed by the collision, 
with the rest-mass density increasing exponentially. The merged object 
expands in a jet-like fashion along the z-direction, with the bulk of the 
matter being accelerated up to 7 ~ 16, or equivalently, v/c 0.998, but 
then settling on much slower fiows with 7 < 2.1. Furthermore, the front of 
the jet has —uq > 1 indicating that part of the shocked matter has sufficient 
energy to have become gravitationally unbound. As a result, the rest-mass 
density at the center of the merged object is smaller than the maximum 
density of the initial configuration, although the origin still represents the 
region where the density is the largest. The right panel, on the other 
hand, refers to a highly-boosted binary, i.e., with vi^/c = 0.8, with each 
star being initially highly distorted by the Lorentz contraction. Also in this 




Figure 4. Representative snapshots of the rest-mass density, p in units where 
c = 1 = Mq (top row), of the Lorentz factor, 7 (middle row), and of the local fluid 
energy, —uq (bottom row), for subcritical models with an initial small boost v\^/c = 0.3 
(left panel) or a large one v\y/c = 0.8 (right panel). Note that the post-collision flow is 
essentially jet-like for the low-boost case (left panel), while essentially spherical for the 
high-boost case (right panel); in this latter case, most of the matter is unbound. 



case, the stars are strongly compressed by the cohision, but the merged 
object expands in a spherical blast-wave fashion, with an almost spherical 
distribution of matter and bulk Lorentz factor. The latter reaches values 
as large as 7 ~ 30, or equivalently, v /c ^ 0.999, which, in contrast with the 
low-boost case, do not decrease in time. As a comparison, the typical bulk 
Lorentz factors obtained in the merger of binary neutron stars in quasi 
circular orbits is 7 1.03 [121]. The very large kinetic energies involved 
in the collision are sufficient to make a very large portion of the stellar 
matter unbound, as clearly shown by the bottom-right panel of Fig. [4j 
which reports the local fluid energy. The rest-mass density distribution 
in the expanding blast wave has a minimum at the origin, where a large 
rarefaction is produced by the matter expanding as an ultrarelativistic 
thick shell. 

The marked transition from a jet-like outflow, not too dissimilar from 
the simple Bjorken flow used to model the very early states of relativistic 
ion-collisions [30 , to a shell-like structure, not too dissimilar from 
"transverse expansion" modelled in the subsequent stages of relativistic 
ion-collisions (see [71j and references therein), signals that it is not 
unreasonable to extrapolate some of the results presented here also to 



the collision of ultrarelativistic elementary particles. 

The transition from the two qualitatively-different regimes discussed 
above is further confirmed by the evolution of the rest-mass normalized 
to the initial value Mq. The simulations in fact reveal that the unbound 
fraction is just a few percent of the total rest-mass in the case of a low-boost 
collision, with most of the matter being confined in the selfgravitating 
"star". This is to be contrasted with what happens for a high-boost 
collision, where the unbound fraction is ~ 100% of the total rest-mass. 
This behaviour provides a strong indication that, at least for subcritical 
collisions, the role played by gravitational forces is a minor one as the 
kinetic energy is increased. This is what happens in the collision of two 
particles at ultrarelativistic speeds, where all of the matter is obviously 
unbound. 

3.3. Critical behaviour and scaling 

A remarkable property of the head-on collision of compact stars is the 
existence of type-I critical behaviour, which was first pointed out in [74] 
and subsequently extended in |75j. In essence, in these works it was 
found that when considering stars with initial zero velocity at infinity, 
it is possible to fine-tune the initial central density (and hence the 
mass) near a critical value, pj, so that stars with pc > p^ would collapse 
eventually to a black hole, while the models with pc < p^ would eventually 
lead to a stable stellar configuration. As a result, the head-on collision of 
two neutron stars near the critical threshold can be seen as a transition in 
the space of configurations from an initial stable solution over to a critical 
metastable one which can either migrate to a stable solution or collapse 
to a black hole [119 . As the critical limit is approached, the survival 
time of the metastable object, Teq, increases as Teq = —A In \pc — with 
A - 10 [111 [75]. 

Although the free-fall velocities considered in [JH [75] were very small, 
the critical behaviour continues to hold also when the stars are boosted to 
ultrarelativistic velocities. Interestingly, the threshold becomes now a 
function of the initial effective boost. Determining becomes especially 
challenging as the Lorentz factor is increased and the dynamics of the 
matter becomes extremely violent, with very strong shocks and rarefaction 
waves. However, it was possible to determine the threshold for all the range 
of initial boosts considered, i.e., V]j/c G [0, 0.95], 75 G [1, 3.2], and even to 
a reasonable accuracy, e.g., = (3.288023 ± 0.000003) x lO^'^g/cm^, for 
the initial boost of v\)/c = 0.3. 

The existence of critical behaviour near which the details of the initial 
conditions become irrelevant and which is the same at different boosts, 
i.e., A does not depend on 7 nor on pc (Refs. [74] 1146] have shown there 



is "universality" when varying 7 and fixing pc), gives us a wonderful tool 
to explore the conditions of black-hole formation also far away from the 
masses and Lorentz factors considered in this paper. This is illustrated in 
Fig. |5j which reports the gravitational mass of the isolated spherical star 
as a function of the effective initial Lorentz factor 

{fdVT^,n^n'^)o' ^ ' 

where T^^y is the stress-energy tensor, is the unit normal to the spatial 
hyperspace with proper volume element dV, and the index refers to 
quantities measured in the initial unboosted frame. I should stress that 
the definition of the effective Lorentz factor (|3| is necessary because the 
stars are extended and thus the Lorentz factor will be different in different 
parts of the star. Expression ([3|, on the other hand, can be taken as ratio 
of the energies measured in the boosted and unboosted frames, and hence 
a generalisation of the Lorentz factor for a point particle (Indeed (7) 1 
for vi) ^ 1). Of course, other parametrizations are possible, still leading to 
scaling laws, but with slightly different exponents. Filled circles indicate 
initial data leading to a black hole, while triangles indicate initial data 
leading to a "star", whereby I mean an object which is at least in part 
selfgravitating (orange errorbars provide an approximate upper limit of 
~ 8% to the error in the measurements). Also indicated as a blue solid 
line is the critical line separating the two regions of black hole and star 
formation (the latter is shown as a shaded region). Clearly, the numerical 
results provide a tight fit of the critical line with a power law 

= ^ ^ 0-92 • (4) 

Mo (7)- (7)1-03 ^ ^ 

Expression Q offers itself to a number of considerations. First, it 
essentially expresses the conservation of energy. Second, in the limit of 
zero initial velocities, (7) ^ 1, one obtains that Mc ^0.92 Mq, so that the 
corresponding total mass, 2Mc, is only ~ 12% larger than the maximum 
mass of the relative spherical-star sequence, i.e., M^g,^ = 1.637 Mq. Third, 
in the opposite limit of (7) ^ 00, expression Q predicts that the critical 
mass will go zero. This is indeed what one would expect: as the kinetic 
energy diverges, no room is left for selfgravitating matter, which will all 
be ejected but for an infinitesimal amount which will go into building the 
zero-mass critical black hole. Fourthly, Q is also in agreement with the 
results in [39, 52 , whereby one can recognize the black- hole formation as 
the crossing of the critical line when moving to larger Lorentz factors while 
keeping the rest-mass constant. 

Finally, using Q it is possible to probe whether the kinetic energies 
achieved by modern particle accelerators, such as the LHC, are sufficient 




to produce micro black holes from the collision of two ultrarelativistic 
particles. Using the results reported in Ref. [69 , the expected energies 
achieved by LHC in the next couple of years will be 4 — 7 TeV, so that a 
proton, whose mass is ~ 938 MeV ~ 8.41 x 10~^^ Mq, can be accelerated 
up to 7 ^ 7.5 X 10^. I have therefore reported the range of masses and 
Lorentz factors accessible to LHC in the inset in Fig. [Sj where it appears 
as a small magenta box. Note that the calculations reported here do not 
intend to be a realistic approximation of the dynamics of ultrarelativistic 
particle collisions. However, these calculations and the presence of a 
critical behaviour can be used to deduce that the ranges reachable at 
the LHC are well below the critical line and thus in the region where a 
partially-confined collided object is expected. 

Of course, this line of arguments wildly extrapolates our results by 
almost 60 orders of magnitude in mass (11 in Lorentz factor) and neglects 
quantum effects and extra-dimension effects that might be important at 
Planck-energy scales. Bearing in mind these caveats, our calculations 
suggest that the production of micro black holes at LHC will be unlikely. 
An additional confirmation that our estimates are not unreasonable comes 



from considering the corresponding energy and Lorentz factors for the 
observed ultra-high energy cosmic rays (UHECR), that are observed with 
energies up to ~ 10^^ eV {i.e., 7 ~ 10^^) and for which there is no evidence 
of black- hole formation when interacting with the atmosphere [32] . Also 
in this case, the relevant range of masses and Lorentz factors is shown in 
the inset and falls in the region where no black holes should be produced. 

As a final remark I note that the scaling relation Q can be expressed 
equivalently in terms of the original stellar compactness, M/R as 

(^] = K' « 0.08 . (5) 

Since Miab = (7)^ is the mass in the lab frame, and since R is the largest 
dimension in that frame being the transverse one to the motion, the ratio 

C^) =K'^^K'^, (6) 
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provides the condition for the amount of energy that, when confined in a 
hoop of radius i?, would lead to a black hole. Hence, expression (|6| extends 
the spirit of the hoop conjecture to the case in which a kinetic energy is 
present. Note that the limiting value (7) = 1 does not corresponds to a 
static configuration (as in the hoop conjecture) but to a binary that is at 
rest at infinity. This explains why in this limit (Miab/^)c = {M/R)c — 
0.08, which is considerably smaller than the value 1/2 predicted by the 
hoop conjecture. 



3.4- Summary 

The calculations reported above demonstrate that it is possible to find 
a criterion for the conditions leading to black-hole formation in the 
collision of two selfgravitating fiuids moving at ultrarelativistic velocities. 
The Lorentz factors reached in these simulations are considerably larger 
than those encountered in merging neutron-star binaries, especially if the 
inspiral is along quasi-circular orbits. The properties of the fiow after the 
collision change with Lorentz factor, with most of the matter being ejected 
in a spherical blast wave for large boosts. Interestingly, the collided object 
exhibits a critical behaviour of type I, which is found to persist also as 
the initial boost is increased. This allows one to derive a simple scaling 
law and extrapolate these results to the energies of elementary particles at 
LHC and conclude that black-hole production is unlikely in that case. 



4. Third piece: horizons as probes of black-hole dynamics 

The third and last "little piece" of numerical relativity that I will 
discuss is instead about calculations in vacuum spacetimes and focuses 



on the merger of two black holes. This process, which represents 
one of the most important source of gravitational waves, is generally 
accompanied by the recoil of the final black hole as a result of anisotropic 
gravitational wave emission. While this scenario has been investigated 
for decades [106 and first estimates have been made using approximated 
and semi-analytical methods such as a particle approximation [55 , post- 
Newtonian methods |147| and the close-limit approximation (CLA) [10], it 
is only thanks to the recent progress in numerical relativity that accurate 
values for the recoil velocity have been computed [24 , 62 , 37, 68 , 80, \S7\ 
[TT01I671I88]. 

Besides being a genuine nonlinear effect of general relativity, the 
generation of a large recoil velocity during the merger of two black holes has 
a direct impact in astrophysics. Depending on its size and its variation with 
the mass ratio and spin, in fact, it can play an important role in the growth 
of supermassive black holes via mergers of galaxies and on the number of 
galaxies containing black holes. Numerical-relativity simulations of black 
holes inspiralling on quasi-circular orbits have already revealed many of 
the most important features of this process showing, for instance, that 
asymmetries in the mass can lead to recoil velocities V]^ < 175 km/s [24[ [62], 
while asymmetries in the spins can lead respectively to ^ 450 km/s or 

^ 4000 km/s if the spins are aligned [68l |80l IllOj or perpendicular to 
the orbital angular momentum [36l [6TJ [37] (see |120j for a review and [89] 
for the most recent results). 

At the same time, however, there are a number of aspects of the 
nonlinear processes leading to the recoil that are far from being clarified 
even though interesting work has been recently carried out to investigate 
such aspects [132 , 93, 82 . One of these features, and possibly the most 
puzzling one, is the generic presence of an ''anti-kick^\ namely, of one (or 
more) decelerations experienced by the recoiling black hole. Such anti- 
kicks take place after a single apparent horizon has been found and have 
been reported in essentially all of the mergers simulated so far (see Fig. 8 
of Ref. [110 for some examples). 

What follows discusses a phenomenological framework which provides 
a novel description of the stages during which the anti-kick is generated, 
and that can be used to formulate a simple and qualitative interpretation 
of the physics underlying this process. I will focus on the head-on collision 
of two nonspinning black holes with different mass. Although this is the 
simplest scenario for a black-hole merger, it contains all the important 
aspects that can be encountered in more generic conditions^ 



Much of what follows is taken from the discussion presented in Refs. [1231 1721 173) . 



4-1- The basic picture 



I will start by presenting a qualitative interpretation of the antikick by 
considering the simple head-on collision of two Schwarzschild black holes 
with unequal masses. This qualitative picture will be made quantitative 
and gauge- invariant by studying the logical equivalent of this process in the 
evolution of a Robinson-Trautman spacetime, with measurements of the 
recoil made at future null infinity. The insight gained with this spacetime 
will be valuable to explain the anti-kick under generic conditions and to 
contribute to the understanding of nonlinear black-hole physics. 

Figure [6] illustrates the dynamics of the head-on collision using a 
schematic cartoon where I am considering a coordinate system centred 
in the total centre of mass of the system and where the smaller black hole 
is initially on the positive z-axis, while the larger one is on the negative 
axis. As the two black holes free-fall towards each other, the smaller 
one will move faster and will be more efficient in "forward-beaming" its 
gravitational wave emission [147j. As a result, the linear momentum will 
be radiated mostly downwards, thus leading to an upwards recoil of the 
black hole binary [cf., stage (1) in Fig. [g]. At the merger, the black-hole 
velocities will be the largest and so will also be the anisotropic gravitational 
wave emission and the corresponding recoil of the system. However, when 
a single apparent horizon is formed comprising the two black holes, the 
curvature distribution on this 2-surface will be highly anisotropic^ being 
higher in the upper hemisphere {cf., red-blue shading in stage (2) of 
Fig. |6|. Because the newly formed black hole will want to radiate all 
of its deviations away from the final Schwarzschild configuration, it will 
do so more effectively there where the curvature is larger, thus with a 
stronger emission of gravitational waves from the northern hemisphere. 
As a result, after the merger the linear momentum will be emitted mostly 
upwards and this sudden change in sign will lead to the anti-kick. The 
anisotropic gravitational wave emission will decay exponentially as the 
curvature gradients are erased and the quiescent black hole reaches its 
final and decelerated recoil velocity [cf., stage (3)Q 

Although this picture refers to a head-on collision, it is supported by 
the findings in the CLA (where the direction of the ringdown kick is 
approximately opposite to that of the accumulated inspiral plus plunge 
kick) [82] and it can be generalized to a situation in which the black holes 

I should remark that other explanations have also been suggested. One of them 
makes use of the Landau-Lifshitz pseudotensor and explains the recoil in terms of the 
cancellation of large and opposite fluxes of momentum, part of which are "swallowed" 
by the black hole 90 . Another one is even more essential and explains the antikick is in 
terms of the spectral features of the signal at large distances, quite independently of the 
presence of a black-hole horizon I113J. All of these views serve the scope of providing 
an intuitive description and are in my view equally valid and useful. 





; .P'CM" 




.R'CM" 










^ P 

few 




(1) 


1 (2) 


(3) 



Figure 6. Cartoon of the generation of the anti-kick in the head-on colhsion of 
two unequal-mass Schwarzschild black holes. Initially the smaller black hole moves 
faster and linear momentum is radiated mostly downwards, thus leading to an upwards 
recoil of the system [stage (1)]. At the merger the curvature is higher in the upper 
hemisphere of the distorted black hole (cf., red-blue shading) and linear momentum is 
radiated mostly upwards leading to the anti-kick [stage (2)]. The black hole decelerates 
till a uniform curvature is restored on the horizon [stage (3)]. 



have different masses, different spins and are merging through an inspiral. 
Also in a more generic case, the newly- formed apparent horizon will 
have a complicated but globally anisotropic distribution of the curvature, 
determining the direction (which is in general varying in time) along which 
the gravitational waves will be emitted. Hence, the geometric properties 
in a dynamical horizon (of a black or white hole) determine its global 
dynamics. I next use the Robinson- Trautman spacetime to validate this 
picture. 

4-2. A useful playground 

The Robinson- Trautman spacetime represents a class of vacuum solutions 
admitting a congruence of null geodesies which are twist and shear- 
free p^6 , with a future stationary horizon and a dynamical past 
(outer trapping) horizon |1Q5[ I143[ 140) l98l llQ8j (past apparent horizon 
hereafter). A Robinson- Trautman spacetime is thus regarded as an isolated 
nonspherical white hole emitting gravitational waves, where the evolution 
of the apparent horizon curvature-anisotropies and the total spacetime 
momentum dynamics can be related unambiguously. The metric is given 



by IS] 



ds' = -(K- ^ - du' - 2dudr + ^dn^ , (7) 

\ r Q J 

where Q = Q{u^Q)^ u is the standard nuh coordinate, r is the affine 
parameter of the outgoing nuh geodesies, and Vt = {6>, 0} are the angular 
coordinates on the unit sphere S'^. Here M^o is a constant and is related 
to the asymptotic mass, while the function K{u^ Q) is the Gaussian 
curvature of the surface corresponding to r = 1 and u = constant, 
K{u, Q) = (5^(1+V^ln(5), where is the Laplacian on S'^. The Einstein 
equations then lead to the evolution equation 

duQ{u, Q) = -Q^\/lK{u, n)/{12Moo) . (8) 

Any regular initial data Q = Q{O^Q) will smoothly evolve according 
to (|8| until it achieves a stationary configuration corresponding to a 
Schwarzschild black hole at rest or moving with a constant speed [41]. 
Equation ([8| implies the existence of the constant of motion A = 

dO./Q'^^ which clearly represents the area of the surface u = const., r = 
const, and can be used to normalise Q so that A = 4:7t. All the physically 
relevant information is contained in the function Q{u, Q)^ and this includes 
the gravitational radiation, which can be extracted by relating Q{u^ Q) to 
the radiative part of the Riemann tensor |431 |12j . 

The past apparent horizon radius R{u, ft) is given by the vanishing 
expansion of the future ingoing null geodesies, satisfying |1Q5| I143j 

Q^Wl \nR = K- 2Moo/R • (9) 

The mass and momentum of the black hole are computed at future null 
infinity using the Bondi four-momentum [91 

with {t]^) = (1, sin^cos(/), sin^sin^, cos6>). Given smooth initial data, 
the spacetime will evolve to a stationary non-radiative solution which, 
in axisymmetry, has the form Q{oo,6) = (1 ^ vx)/Vl — v'^, with 
X = cosO [91]. The Bondi four-momentum associated to Q{oo^6) has 
components 

(P(oo))" = (Moo/\/1 - «2) (1, 0, 0, ±v) , (11) 

so that the parameter v in Q{oo^6) can be interpreted as the velocity of 
the Schwarzschild black hole in the z-direction. 



One of the difficulties with Robinson- Trautman spacetimes is the 
definition of physicahy meaningful initial data. Although this is meant 
more as a proof-of-principle than a realistic configuration, it is possible to 
adopt the prescription suggested in Ref. [ 12| 



g(o,^) = Qo 



y^^Y^^UJX Vl + WX_ 



(12) 



and which was interpreted to represent the final stages {i.e., after a 
common apparent horizon is formed) of a head-on collision of two boosted 
black holes with opposite velocities w and mass ratio q [12 . In practice, to 
reproduce the situation shown in Fig. [6j it is sufficient to choose w < and 
take q G [0^ 1]; a more general class of initial data and the corresponding 
phenomenology can be found in [72l [13]. Note that Qo is chosen so 
that to ^ = 47r and that in general the deformed black hole will not be 
initially at rest. As a result, given the initial velocity vq = P^(0)/P^(0), 

a boost is performed transformation = A"^(i'o)P^ so that P (0) = 
by construction. The numerical solution of eq. ([8| with initial data ( 12 ) is 
performed using a Galerkin decomposition as discussed in detail in [91 . 

Figure [7| reports the typical evolution of a Robinson- Trautman 
spacetime with the lower panel showing the evolution of the curvature 
of the past apparent horizon Kj^^ = 2M^/ R^{x) at the north {x = 1) and 
south pole {x = —1), and with the upper panel showing the evolution of the 
recoil velocity. Note that the two local curvatures are different initially, 
with the one in the upper hemisphere being larger than the one in the 
lower hemisphere {cf., Fig. [6|. However, as the gravitational radiation 
is emitted, this difference is erased. When this happens, the deceleration 
stops and the black hole attains its asymptotic recoil velocity. The inset 
reports the curvature difference relative to the asymptotic Schwarzschild 
one, Kj^^ — 1, whose exponentially decaying behaviour is the one expected 
in a ringing black hole (see also Fig. 7 of Ref. [72]). 

As mentioned before, the one shown in Fig. [7] is a typical evolution of a 



Robinson- Trautman spacetime and is not specific of the initial data (12). 
By varying the values of k;, in fact, it is possible to increase or decreases 
the final recoil and a sign change in w simply inverts the curvature at the 
poles so that, for instance, initial data with w > would yield a black hole 
accelerating in the positive z-direction. Interestingly, it is even possible 
to fine-tune the parameter w so that the recoil produced for a Robinson- 
Trautman spacetime mimics the anti-kick produced by the quasi-circular 
inspiral of nonspinning binaries. This is shown in Fig. [T] which reports 
the recoil as a function of the symmetric mass ratio z^ = g/(l + g')^, and 
where the dashed line refers to the anti-kick for the inspiral of nonspinning 
binaries in the CLA [82] (the parameters chosen, i.e., w = —0.425 and 
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Figure 7. Typical evolution of a Robinson- Trautman spacetime. Shown in the 
lower panel is the evolution of the curvature i^^H north (x = 1) and south pole 

(x = —1). Shown in the upper panel is the evolution of the recoil, which stops decreasing 
when the curvature difference is erased by the emitted radiation (dotted line). Note 
that the curvature decays exponentially to that of a Schwarzschild black hole (inset). 

ri2 = 2 M, are those minimising the differences). Considering that the 
two curves are related only logically and that the CLA one contains all the 
information about inspiralling black holes, including the orbital rotation, 
the match is surprisingly good. 

It is also suggestive to think that the curve in Fig. |8] is actually 
composed of two different branches, one of which is characterized by large 
curvature gradients across the apparent horizon but small values of the 
curvature (this is the \ow-u branch and is indicated with squares), while 
the other is characterized by small curvature gradients and large values of 
the curvature (this is the high-i^ branch and is indicated with circles). The 
same recoil velocity can then be produced by two different values of for 
which the effects of large curvature gradients and small local curvatures 
are the same as those produced by small curvature gradients but large 
local curvatures. 

To go from this intuition to a mathematically well-defined measure 
one can compute the mass multipoles of the intrinsic curvature of the 
initial data using the formalism developed in [T5^ for dynamical horizons. 





Figure 8. Recoil velocity shown as a function of the symmetric mass ratio v when 
w = —0.425, with the dashed line refers to the anti-kick from the inspiral of nonspinning 
binaries in the CLA [82]. Note that the curve can be thought as composed of two 
different branches. 



Namely, it is possible to compute the mass moments as (the mass-current 
are obviously zero) 

where Pn{x) is the Legendre polynomial in terms of the coordinate x{0) 
which obeys dex = - sin (9i?(i9)V(i?2^Q((9)2), with R^^ = y^A^^/{A7r) 
and x{0) = 1. Using these multipoles it is possible to construct an effective- 
curvature parameter Kq{^ that represents a measure of the global curvature 
properties of the initial data and from which the recoil depends in an 
injective way. Because this effective-curvature parameter has to contain 
the contribution from the even and odd multipoles, the expression 

i^eff = M2|^M2„+i/3"-i|, (14) 

n=l 



was found to reproduce exactly what expected (note Mi = to machine 
precision). 




Figure 9. Recoil velocity shown as a function of the effective curvature. In contrast 
with Fig. [s] which uses the same symbols employed here, the relation between the 
curvature and the recoil is now injective. 



This is shown in Fig. |9j which reports the recoil velocity as a function 
of iCeff- As predicted, and in contrast with Fig. |8j the relation between 
the curvature and the recoil is now injective, with the maximum recoil 
velocity being given by the maximum value of i^eff (see inset), and with 
the two branches coinciding. The expression ( [m] ) suggested above for 
i^Teff is not unique and indeed a more generic one will have to include 
also the mass-current multipoles to account for the spin contributions (see 
discussion below). However, lacking a rigorous mathematical guidance, 
the phenomenological i^eff is a reasonable, intuitive approximation. 



4.3. A more general view 

Despite the valuable insight, the treatment summarised above and 
presented in Ref. |123j had obvious limitations. First, the Ansatz (14) 
for Kes, i.e., i^eff = /even (^2^) X /odd (^2^+1) IS uot Straightforwardly 
generalize to the non-axisymmetric case. Second, the functions /even and 
/odd can be written in the simplest possible form, i.e., as a linear expansion 
in Mi's, i.e., Keff = (a2M2 + a4M4 + . . .) x {asMs + a^M^ + • • where 
the phenomenological coefficients a^'s depend on the details of the 



employed initial data. Finally, the white-hole horizon analysis in Robinson- 
Trautman spacetimes needs to be extended to the genuine black-hole 
horizon case. 

While the focus in what discussed above (and presented in Ref. |123j ) 
was on expressing the difference between the final kick velocity Voo and 
the instantaneous kick velocity V]^{u) at an (initial) given time in terms 
of the geometry of the common apparent horizon at that time it is 
possible to derive a more generic view based on geometric quantities that 
are evaluated at a given time during the evolution. More specifically, it is 
possible to consider the variation of the Bondi linear momentum vector in 
time {dPf / du){u) as the relevant geometric quantity to monitor at null 
infinity . This quantity can then be correlated with a counterpart on 
the black- hole horizon H^, e.g., a vector Kl^{v) (function of an advanced 
time v), which represents an extension of the effective curvature introduced 
in the previous Sectiorj^ 

In the case of a Robinson- Traut man spacetime, the causal relation 
between the white-hole horizon 1-L~ and null infinity made possible to 
establish an explicit functional relation between dv]^/du and K'^^{u). In the 
case of a generic black-hole horizon, such a direct causal relation between 
the inner horizon and is lost. However, since the corresponding causal 
pasts of ^+ and 1-L~ coincide in part, non-trivial correlations are still 
possible and expected. These correlations can be measured by comparing 
geometric quantities /imn('^) at and ho^tiu) at J^~^^ both considered 
here as two timeseries^ In particular, it is reasonable to take Klf^{v) as 
hinn{v) and {dPf / du){u) as ho^t{u). 

This approach resembles therefore the methodology adopted in 
scattering experiments. Gravitational dynamics in a given spacetime 
region affects the geometry of appropriately-chosen outer and inner 
hypersurfaces of the black-hole spacetime. These hypersurfaces are then 
understood as test screens on which suitable geometric quantities must 
be constructed. The correlations between the two encode geometric 
information about the dynamics in the bulk, providing information 
useful for an inverse- scattering approach to the near- horizon dynamics. 
In asymptotically flat black-hole spacetimes, null infinity and the 
(event) black-hole horizon provide natural choices for the outer and 
inner screens. This is summarised in the Carter-Penrose diagram in 

^ Another appealing approach that has a similar goal of correlating strong-fields effects 
with (the visualization of) spacetime curvature has been proposed recently by the group 
in Caltech p^[99l. 

^ Note that the meaningful definition of timeseries cross-correlations requires the 
introduction of a (gauge-dependent) relation between advanced and retarded time 
coordinates v and u. In an initial value problem this is naturally provided by the 
3+1 spacetime slicing by time t. 



Figure 10. Carter-Penrose diagram illustrating the scattering approach to near- 
horizon gravitational dynamics in a generic spherically symmetric collapse. The 
event horizon Ti'^ and null infinity provide spacetime canonical screens on 

which geometric quantities, respectively accounting for horizon deformations and wave 
emission, are defined. Their cross-correlation encodes nontrivially information about 
the bulk spacetime dynamics. 



Fig. 10, which illustrates the cross-correlation approach to near-horizon 
gravitational dynamics. The event horizon and null infinity J^~^ 
provide spacetime screens on which geometric quantities, accounting 
respectively for horizon deformations and wave emission, are measured. 
Their cross-correlation encodes information about the bulk spacetime 
dynamics. 

The picture offered by Fig. [lO] can be easily adapted to the 3+1 
approach commonly adopted in numerical relativity. Since neither the 
black-hole event horizon nor null infinity are in general available during the 
evolutiorp^ it is possible to adopt as inner and outer screens a dynamical 
horizon (future outer trapping horizon [66l [161 Hi]) timelike 
tube B at large spatial distances, respectively. In this case, the time 
function t associated with the 3+1 spacetime slicing provides a (gauge) 



The latter would properly require either characteristic or a hyperboloidal evolution 
approach. 



r = 








Figure 11. Carter-Penrose diagram for the scattering picture in a Cauchy initial 
value approach. The dynamical horizon 7^+ and a large-distance timelike hypersurface 
B provide inner and outer screens. Note that the dynamical horizon is split in two 
portions: outer and inner (solid and dashed blue lines, respectively) and that the 3+1 
slicing sets a common time t for cross-correlations. 



mapping between the retarded and advanced times u and so that cross- 
correlations between geometric quantities at and B can be calculated as 
standard timeseries hinn{t) and /lout(^)- This is summarised in the Carter- 
Penrose diagram in Fig. [TT] which is the same as in Fig. 10, but where the 
3+1 slicing sets an in-built common time t for cross-correlations between 
the dynamical horizon {i.e., the inner screen) and a large-distance 
timelike hypersurface B {i.e., the outer screen). 

Within this conceptual framework it is then possible to define a 
phenomenological curvature vector Kf^{t) in terms of the mass multipoles 
of the Ricci scalar curvature at and show that this is closely 
correlated with a geometric quantities {dP^ /dt){t)^ representing the 
variation of the Bondi linear momentum time on . How to do this 
in practice for a black-hole spacetime requires much more space that I 
can take in this contribution and therefore refer the interested reader to 
Refs. [Ill [73], where this is discussed in great detail. 



4-4' Summary 

The discussion reported above demonstrates that quahtative aspects of 
the post-merger recoil dynamics at infinity can be understood in terms of 
the evolution of the geometry of the common horizon of the resulting 
black hole. Moreover, suitably-built quantities defined on inner and 
outer worldtubes (represented either by dynamical horizons or by timelike 
boundaries) can act as test screens responding to the spacetime geometry 
in the bulk, thus opening the way to a cross-correlation approach to probe 
the dynamics of spacetime. This picture was shown to hold both for a 
simple Robinson- Trautman spacetime, but also for more generic binary 
black-hole spacetimes. In this latter case, this is possible through the 
construction of a phenomenological vector Kf^{t) from the Ricci curvature 
scalar on the dynamical horizon sections, which then captures the global 
properties of the fiux of Bondi linear momentum [dPf / dt){t) at infinity, 
namely the acceleration of the BH. 

A geometric framework looking at the horizon's properties offers a 
number of connections with the literature developing around the use 
of horizons to study the dynamics of black holes, as well as with the 
interpretations of such dynamics in terms of a viscous-hydrodynamics 
analogy. Much of the machinery developed using dynamical trapping 
horizons as inner screens can be extended also when a common horizon is 
not formed (as in the calculations reported in Ref. |139j ). While in such 
cases the identification of an appropriate hypersurface for the inner screen 
can be more difficult, once this is found its geometrical properties can be 
used along the lines of the cross-correlation approach discussed here for 
dynamical horizons. 

5. Conclusions 

The "three little piece" for numerical computer and relativity presented in 
the sections above ought to provide a reasonable idea of the "Renaissance" 
that numerical relativity is now experiencing. More importantly, they 
should be able to convey the enormous potential that numerical-relativity 
simulations have in revealing aspects of the theory that cannot be handled 
analytically, or in exploring nonlinear regimes that cannot be investigated 
through perturbative approaches. As remarked repeatedly, the examples 
brought represent only a personal (and biased) selection of the intense 
work carried out recently and surely are not exhaustive in terms of the 
physical scenarios that can be explored. Much more can be said about 
this and surely it will not have to wait for the bicentenary of Einstein's 
stay in Prague. 
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